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Abstract. We have investigated a group of stars known to have low chromospheric ages, but high kinematical ages. 
Isochrone, chemical and lithium ages are estimated for them. The majority of stars in this group show lithium 
abundances much smaller than expected for their chromospheric ages, which is interpreted as an indication of 
their old age. Radial velocity measurements in the literature also show that they are not close binaries. The results 
suggest that they can be formed from the coalescence of short-period binaries. Goalescence rates, calculated taking 
into account several observational data and a maximum theoretical time scale for contact, in a short-period pair, 
predict a number of coalesced stars similar to what we have found in the solar neighbourhood. 
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1. Introduction 

The chromospheric activity of a late-type star is frequently 
interpreted as a sign of its youth. Young dwarfs show high 
rotation rates, and the interaction between rotation and 
outer envelope convection is expected to drive the chromo¬ 
spheric activity. Nevertheless, not only young single stars 
present high rotation rates. Close and contact binaries can 
keep high rotation over several billion years. In such stars, 
the rotational angular momentum loss is balanced by the 
proximity of the stars in the system, which results in trans¬ 
fer of orbital angular momentum to the rotational spins. 
Thus in chromospheric activity surveys aimed at late type 
stars, we expect to find two classes of objects: young stars 
and chromospherically active binaries. Sometimes, a star 
suspected of being young can b e instead a spectroscopic 
binary (Soderblom et al. 1998|) , not yet investigated by 
radial velocity surveys. 

The chromospheric activity surveys by the Mount 
Wilson group are directed at, but not only at, late type 
stars. The two surveys that comprise the bulk of a sample 
used by some of us in the derivation of chemodynamical 
constraints on the evolution of the Galaxy (Rocha-Pinto 
et al. ^OOOa , 2000b ) were based on solar-type stars, in a 
spectral range from F8 V to K4 V. Chromospherically 
active binaries were generally avoided, since the surveys 
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investigate the chromospheric activity in single stars. Due 
to this, the division of these surveys into two classes, of 
active and inactive stars, corresponds closely to an age 
segregation. 

This was very well demonstrated by Soderblom ( |1990| , 
see also Jeffries & Jewell 1993), who studied the kinemat¬ 
ics of active and inactive stars. The active stars are con¬ 
centrated in a region of low velocities in a space velocity 
diagram, as expected for young objects; on the other hand, 
the inactive stars are scattered in this diagram, just like 
old stellar populations. Few active stars do not follow this 
rule, showing considerably high velocities. Soderblom calls 
attention to them, but int erpret them as possible runaway 
stars. Rocha-Pinto et al. ( 2002 ) have increased the num¬ 
ber of active stars with spatial velocities to 145. Several 
of these stars, show velocities which are inconsistent with 
their presumed age. 

The term CYKOS (acronym for chromospherically 
young, kinematically old stars) is applied here to all chro¬ 
mospherically active stars which, in a velocity diagram, 
present velocity components greater than the expected 
value for such stars, irrespective of the fact that this object 
is an undiscovered close binary, a runaway star or another 
kind of object. 

This paper presents several newly identified CYKOS, 
and proposes an explanation for some of them. It is orga¬ 
nized as follows: in section 2, we present the sample and 
the method used to define a CYKOS. Section 3 analy- 
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ses what is presently known about their ages from chro¬ 
mospheric and isochrone age measurement methods. In 
section 4, we show that these objects have very low Li 
abundances compared to other stars with the same tem¬ 
perature. A critical review of the literature about some 
individual CYKOS follow in section 5. Finally, in section 
6, we propose that a small number of these objects could 
probably have been formed by coalescence of binaries. 


2. Identification of CYKOS 


The CYKOS can be identified by a diagram of spatial ve¬ 
locities {U xV otW xV) showing only active stars, in 
analogy to their first discovery by Soderblom. Not all ob¬ 
jects identified 'm aU xV diagram are also identified in a 
W X V diagram, and vice versa. We expect that CYKOS 
showing high velocities in more than one component are 
really peculiar objects, and not just stars having a com¬ 
ponent velocity in the tail of the distribution. 

Some CYKOS also appear clearly in an age-velocity 
diagram, so that, in principle, this diagram could also be 
used to identify them. However, the uncertainty in the 
chromospheric ages can make some inactive stars appear 
young in such a plot, and our main purpose is not the 
study of young inactive stars, but rather of active stars 
that are kinematically old. 

We have used as our primary source the samp le con - 
taining 145 active stars built by Rocha-Pinto et al. ( 20021 ). 
This sample is composed by all active stars from Rocha- 
Pinto et al. ( 2000bD for which radial velocity measure¬ 
ments are available in the literature. The heliocentric 
velocities are calculate d with the equations provided by 
Johnson & Soderblom ( 1987 ). 

Figs.j^ and b show the U x V and W x V diagrams 
for these stars. The velocities were corrected for the so¬ 
lar motion (17 q,Vo,1F q) = (—9,11,6) km s“^ according 
to Mihalas & Binney ( 198l| ). In the plots, the semiaxes 
of the inner and outer ellipses are equal to la and 3cr, 
respectively, where a is the velocity dispersion of stars 
with ages lower than 1 Gyr. We have adopted au = 20 
km/s, av = 10 km/s and aw — 8 km/s, which corre¬ 
spond to the velocity dispersion of the youngest stellar 
popul ation (ages 0 to 1 Gyr), according to Rocha-Pinto 
et al. ( 2002| ), and have considered as CYKOS all objects 
located beyond the outer ellipses. The velocity dispersions 
for the youngest population are in good agreement with 
other determinations in the literature. The values for au 


range from 11.7 km/s (Meusinger et al. 1991) to 23 km/s 
(Cayrel de Strobel 1974). The agreement is closer for ay, 
ranging from 5.9 km/s (Meusinger et al. 1991) to 10.1 
km/s (Cayrel de Strobel 1974), and ctw, ranging from 8 
km/s (Wielen |1974|) to 9.7 km/s (Cayrel de Strobel |1974|). 

In Tab. (0)7 we list the thirty stars identified as 
CYKOS. The first column gives the name of the star, 
followed by the spectral type, [B — V) colour, chromo¬ 
spheric activity index, chromospheric and isochrone ages, 
as commented upon later, and the U, V and W heliocen¬ 
tric velocities, which are given with formal uncertainties. 


The last column shows the velocity criteria used in order 
to classify the object as a CYKOS. 

The results also show that CYKOS are generally more 
apparent in V than in the other velocity components, as 
expected for a kinematically old stellar population: the 
V component is systematically more negative than in the 
case of normal active stars, which are presumably young. 
This is an effect of the asymmetrical drift. The stars are 
likely to acquire increasing random velocities with respect 
to the local standard of rest, due to subsequent encounters 
with giant molecular clouds. In U and W, there will be 
a symmetric increase in the velocity dispersion and we 
would not expect the CYKOS to be much different from 
the normal stars, from the consideration of these velocity 
components only. 

The peculiar character of these objects can be seen 
from Figs. ^ and d, where the same diagrams are shown 
for 102 3 A d warfs, taken from the large compilation by 
Palous ( |1983| ). Since A dwarfs are very young stars, due 
to their maximum life expectancy, their kinematical prop¬ 
erties are consistent with those of very young late-type 
dwarfs. We do find stars having velocity components 
greater than 3a, where we have used the same velocity 
dispersions used in Figs, pi and b. Nevertheless, these A 
dwarf outliers have smaller velocities than the average ve¬ 
locities of the CYKOS as shown in Figs. |^a and b. Aroung 
4% of the A dwarfs are outliers, while this number goes to 
20% in the case of the late-type dwarfs. For the A stars, 
these outliers probably represent the tail of the velocity 
distribution, which is reinforced by the fact that their V 
distribution is nearly symmetrical, contrary to what hap¬ 
pens for the CYKOS. 


3. What age measurement methods tell us 

The chromospheric ages of half of the objects listed in 
Tab. (|^) are lower or similar to 2 Gyr. This can be seen 
from the fourth and fifth columns of the Table, which give 
the chromospheric activity index log i?HK chro¬ 

mospheric ages in Gyr, respectively. The chromospheric 
index was defined by Noyes et al. ( 1984 ) and the calcula¬ 
tion of chr omos pheric ages was discussed by Rocha-Pinto 
& Maciel ( 1998 ). 

There are 9 stars in Table ^ whose chromospheric age 
is larger than 3 Gyr. Strictly speaking, they cannot be 
considered as ‘chromospherically young’, and are listed in 
view of their activity levels (logi?(jj^ > —4.75), which tra¬ 
ditionally indicate a young age (Soderblom 199C| ). Five 
of them have logi?(jj^ < —4.70. Given that the error in 
logi?(jj^ is expected to be around 0.04 dex, it is possible 
that these stars are inactive, rather than active stars. On 
the other hand, at least one of them presents significant X- 
ray emission (HD 89995). Also, we must take into account 
that a high chromospheric age could be caused by the 
metallicity of the star. Metal-poor CYKOS are also chro¬ 
mospherically older, due to the metallicity depe ndenc e of 
the chromospheric age (Rocha-Pinto & Maciel 1998 ). In 
fact, some of these 9 stars are metal-poor, in compari- 
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Table 1. Identified CYKOS. The columns list the stellar name, MK type, {B — V), chromospheric index, chromospheric 
and isochrone ages, heliocentric spatial velocities U, V and W, and the components, according to which the object 
was classified as CYKOS. In the sixth cloumn, the remarks ‘ZAMS’, ‘red MS’ and ‘sd’ refer, respectively, to stars near 
the zero age main sequence, in the red part of the main sequence, or that are subdwarfs according to their position, 
and for which no age determination was possible. 


HD/BD 

MK 

(B-V) 

log 

chrom.age 

isoch.age 

U 

V 

W 

criteria 

5303 

G3: V-t 

0.71 

-4.03 

0.38 

2.30 

81 ±5 

-58 ±3 

-13 ± 1 

UV 

7983 

G2 V 

0.59 

-4.75 

7.91 

18.0 

-138 ± 11 

-72 ±6 

-51 ±5 

UVW 

13445 

K1 V 

0.77 

-4.74 

2.34 

redMS 

98 ± 0.7 

-75 ±1 

-25 ±2 

UVW 

16176 

F5 V 

0.48 

-4.73 

2.07 

2.0 

25 ±2 

-46 ±2 

-22 ± 1 

V 

20766 

G2.5 V 

0.64 

-4.65 

2.09 

redMS 

71 ±0.5 

-47 ± 1 

16 ± 1 

UV 

39917 

G8 V 

0.76 

-4.05 

0.28 

1.3 

-61 ±6 

-31 ±2 

-19 ±2 

U 

51754 

GO 

0.57 

-4.56 

3.35 


-194 ± 17 

-139 ± 13 

3± 1 

UV 

65721 

G6 V 

0.74 

-4.67 

2.86 

redMS 

-37 ±2 

-39 ±2 

28 ± 1 

V 

74385 

K1 V 

0.91 

-4.55 

0.93 

redMS 

15 ±0.4 

-21 ±2 

-30 ±0.5 

W 

88742 

G1 V 

0.62 

-4.69 

2.53 

10.0 

36 ± 0.5 

-45 ± 1 

-4± 1 

V 

89995 

F6 V 

0.46 

-4.74 

4.70 

2.6 

55 ±2 

-42 ±2 

-17±3 

V 

103431 

dG7 

0.76 

-4.68 

1.74 

redMS 

68 ±4 

-39 ±2 

-13 ±2 

UV 

106516 

F5 

0.46 

-4.65 

6.43 

7.4 

-54 ± 1 

-74 ±2 

-59 ±2 

VW 

120237 

G3 IV-V 

0.58 

-4.75 

4.16 

13.8 

47 ±2 

-57 ±2 

-3± 1 

V 

123651 

GO/Gl V 

0.53 

-4.74 

5.39 

13.0 

29 ±4 

-10 ±4 

-33 ±2 

W 

131582 

K3 V 

0.96 

-4.73 

2.79 

redMS 

66 ±2 

-67 ±2 

15 ±2 

UV 

131977 

K4 V 

1.11 

-4.49 

0.21 

redMS 

-49 ±2 

-22 ± 1 

-32 ± 1 

W 

144872 

K3 V 

0.96 

-4.74 

2.89 

redMS 

-70 ± 1 

1± 1 

-2 ±2 

u 

149661 

K2 V 

0.81 

-4.58 

0.83 

redMS 

2±2 

-1±0.4 

-30 ± 1 

w 

152391 

G8 V 

0.76 

-4.39 

0.28 

redMS 

-85 ±2 

-112 ±2 

8± 1 

UV 

165401 

GO V 

0.63 

-4.65 

3.92 

18.9 

79 ±2 

-90 ± 1 

-40 ± 1 

UVW 

183216 

G2 V 

0.60 

-4.62 

0.50 

1.8 

40 ±2 

-46 ± 1 

0± 1 

V 

189931 

G1 V 

0.60 

-4.64 

0.76 

ZAMS 

41 ±2 

-50 ± 1 

1± 1 

V 

196850 

GO 

0.57 

-4.64 

1.31 

redMS 

0±0.4 

-22 ±2 

-33 ± 1 

w 

209100 

K4.5 V 

1.06 

-4.56 

0.39 

redMS 

80 ± 1 

-41 ±0.2 

4± 1 

UV 

219709 

G2/G3 V 

0.65 

-4.62 

1.14 

8.5 

25 ± 1 

-42 ± 1 

-10 ±0.1 

V 

230409 

GO 

0.70 

-4.70 

8.78 

sd 

-133 ± 16 

-129 ± 13 

-13 ± 1 

UV 

-tl5 3364 

GO 

0.63 

-4.43 

0.49 

8.7 

-59 ±3 

-22 ±2 

2± 1 

u 

+51 1696 

sdGO 

0.55 

-4.42 

3.51 


228 ± 23 

-277 ± 33 

53 ± 1 

UVW 


son with the majority of the other stars. However, since 
the metallicity introduces an additional source of error in 
the calculation of the chromospheric age, it is not unlikely 
that these objects could be younger than what is shown in 
Table |^. For instance, we have remarked that even some 
of these older CYKOS have velocities considerably larger 
than the mean velocity of their coeval stars (this is par¬ 
ticularly true for the CYKOS having chromospheric age 
between 2 and 4 Gyr). For these reasons only, we have 
decided to keep them in the sample. 

The basic problem deserving an explanation is why the 
chromospheric ages of these objects are low (sometimes 
very low), while their kinematic ages are high (sometimes 
very high). What can be said about their ages from other 
methods? 

The objects in Tab.(|^ have a broad metallicity dis¬ 
tribution. Most of these stars have [Fe/H] between —0.40 
and -1-0.20, but 11% of them have photometric metallici- 
ties lower than —0.60 dex. From the point of view of chem¬ 
ical evolution, they have a broad age range, with averages 


around 3-5 Gyr if we adopt the age-metalliticy relation 
given by Rocha-Pinto et al. ( ^000b| ). 

For the calculation of isochro ne ag es, we have used 
the isochrones by Bertelli et al. (1994). The age of the 
star was calculated by interpolation in each isochrone. A 
final interpolation, taking into account the ages at several 
metallicities (that is, several isochrone grids), uses the real 
stellar metallicity to find the stellar age. The results are 
shown in the sixth column of Tab. 

These ages are somewhat rough and some care must 
be taken when interpreting these results. The reason for 
this is that the mi deficiency, present in active stars (see 
Rocha-Pinto & Maciel 1998, and references therein), hin¬ 
ders the determination of accurate stellar parameters by 
photometric indices. In our case, logTeff, which is esti¬ 
mated from the equations given by Olsen ( 1984 ), can be 
miscalculated. Three stars do not have parallaxes mea¬ 
sured by HIPPARGOS, or have uvby colours outside the 
range covered by the calibrations by Olsen, and do not 
have isochrone ages in Table |^. 
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Fig. 1 • Panels a and b: spatial veloc¬ 
ity diagrams for 145 active stars. The 
ellipses have semiaxes equal to a and 
3a, where a is the velocity dispersion 
for stars with ages lower than 1 Gyr 
(Rocha-Pinto et al. 2002 ) in the cor¬ 
responding velocity components. For 
these dispersions, we use (t„ = 20 km/s, 
Gv = 10 km/s and Uu, = 8 km/s. 
Objects located beyond the outer el¬ 
lipses, corresponding to 3a, are consid¬ 
ered CYKOS. Panels c and d: spatial 
velocity diagrams for 1023 A dwarfs , 
taken from the sample of Palous (1983). 
In all panels, we have used spatial ve¬ 
locities corrected f or the solar motion 
(Mihalas & Binney 1981). 


Nearly 40% of the stars lie in the red part of the main 
sequence, for which no age determination is possible. The 
remaining stars are distributed nearly equally between 
young (7 stars having less or about 2 Gyr) and old stars (8 
stars with ages greater than 7 Gyr). There are no preferred 
ages for these stars. 

Our results show that, in spite of some of CYKOS hav¬ 
ing low isochrone ages, others can be very old. The cooler 
stars in the red MS can have very different ages, since they 
are in a colour range where no perceptible evolution in the 
HR diagram is visible. 

4. Lithium in CYKOS 

Although Li depletion and production in stars are pro¬ 
cesses not completely well understood, in some cases the Li 
abundance could be used as a youth indicator. If CYKOS 
are young objects, as suggested by their chromospheric ac¬ 


tivity, they must present high Li abundances. On the other 
hand, if they show very depleted Li, they must be evolved 
objects, irrespective of what their magnetic activity might 
tell us. 

4.1. Observations 

We have obtained spectra for 28 stars, including CYKOS 
and normal active stars (used as reference objects). The 
observations were carried out at the Laboratorio Nacional 
de Astrofisica (LNA, Brazil) in two observing runs, in 
August 1998 and July 1999. The spectra were obtained 
with the coude spectrograph at the 1.6m telescope, using 
a SITe CCD of 1024x1024 pixels with 24/imx24/j,m pixel 
size, and a grating of 1800 1/mm yielding a resolution over 
two pixels of w 24,000 covering the wavelength range AA 
6640-6780 A. 
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Table 2. Observations of CYKOS and active stars. 


HD/BD 

R.A 

h 

. (2000) 

m s 

Dec. (2000) 

o f n 

J.D. 

(2451000 +) 

exp. 

(s) 

Teff 

(K) 

logff 

(km/s) 

[Fe/H] 

N(Li) 

Vr 

(km/s) 

870 

00 

12 

51 

-57 

54 

48 

037.782 

900 

5350 

4.40 

0.5 

-0.12 

0.50 

-2 

1237 

00 

16 

04 

-79 

51 

02 

384.785 

660 

5250 

4.33 

0.7 

0.00 

1.90 

-3 

13445 

02 

10 

15 

-50 

50 

00 

037.836 

900 

5270 

4.50 

0.8 

-0.10 

0.50 

+53 

17051 

02 

42 

31 

-50 

48 

12 

384.837 

450 

6040 

4.25 

1.0 

0.10 

2.40 

+4 

20766 

03 

17 

36 

-62 

35 

04 

037.850 

600 

5715 

4.30 

1.0 

-0.20 

0.50 

+9 

22049 

03 

32 

59 

-09 

27 

31 

384.832 

220 

5215 

4.83 

1.3 

-0.15 

0.75 


106516 

12 

15 

10 

-10 

17 

54 

384.398 

600 

6190 

4.20 

1.0 

-0.55 

< 1.20 

+8 

124580 

14 

15 

38 

-44 

59 

55 

037.410 

1200 

5845 

4.28 

1.5 

-0.20 

2.75 

+4 

131977 

14 

54 

32 

-21 

11 

28 

033.452 

900 

4585 

4.58 

1.0 

0.20 

< -1.00 

+26 

138268 

14 

15 

58 

-44 

59 

55 

033.502 

2400 

5975 

4.40 

1.0 

0.30 

2.90 

-59 








037.426 

1200 





2.90 


149661 

16 

36 

19 

-02 

19 

13 

037.442 

600 

5235 

4.50 

1.0 

0.05 

< -1.00 

-14 

152391 

16 

53 

01 

-00 

00 

22 

037.455 

900 

5450 

4.37 

1.2 

0.00 

1.10 

+44 

154417 

17 

05 

16 

00 

42 

25 

037.481 

666 

5970 

4.35 

1.3 

-0.10 

2.80 

-18 

165401 

18 

05 

37 

04 

39 

42 

037.468 

900 

5755 

4.21 

1.0 

-0.40 

0.50 

-120 

174429 

18 

53 

05 

-50 

10 

49 

037.528 

1320 

5100 

4.10 



3.00 

-10 

181321 

19 

21 

29 

-34 

58 

56 

037.567 

780 

5975 

4.30 

1.5 

0.10 

3.10 

-19 

185124 

19 

37 

46 

-04 

38 

48 

037.579 

480 

6760 

4.20 



2.80 

-33 

189931 

20 

04 

02 

-37 

52 

15 

037.587 

900 

5865 

4.35 

1.0 

0.10 

2.15 

-44 

202628 

21 

18 

25 

-43 

20 

05 

384.697 

660 

5750 

4.24 

1.0 

0.05 

2.15 

+9 

202917 

21 

20 

49 

-53 

01 

58 

384.672 

1800 

5555 

4.25 

1.6 

0.05 

3.35 

-1 

206667 

21 

44 

44 

-42 

07 

46 

384.709 

1000 

5950 

4.22 

1.0 

0.00 

2.40 

+ 14 

209100 

22 

03 

21 

-56 

47 

09 

384.664 

400 

4660 

4.90 

1.5 

0.10 

0.15 

-44 

217343 

23 

00 

18 

-26 

09 

05 

384.723 

900 

5755 

4.44 

1.3 

0.00 

3.20 

-9 

221231 

23 

31 

00 

-69 

04 

29 

384.751 

800 

5910 

4.42 

1.3 

0.00 

2.95 

+2 

223913 

23 

53 

40 

-65 

56 

55 

384.762 

660 

5985 

4.50 

1.0 

0.15 

2.65 

+18 

+15 3364 

18 

07 

18 

15 

56 

54 

037.491 

1800 

5685 

4.20 

0.8 

0.00 

0.00 

+24 

55 







037.504 

1800 





0.00 



The data were reduced using standard tasks of IRAF 
package. Spectra from different exposures were added by 
weighting them with (S/N)^. The final S/N ratio for the 
stars are in the range of 100 to 200. The spectrum of a hot 
star obtained with the same configuration was inspected 
for telluric lines. The log of observations is reported in 
Table ^ together with some other information. The radial 
velocities listed were calculated using a set of unblended 
atomic lines in the same data. The error in the radial 
velocity is 2 km/s. 

4.2. Spectra of normal stars and CYKOS 

In Fig. we show spectra in the Li region (A 6707 A) 
for some objects taken as normal. The spectra in this plot 
are ordered according to the chromospheric activity, the 
lowest spectra being that of the least active star. 

The chromospheric activity order must represent ap¬ 
proximately an age order: the most active objects are sup¬ 
posed to be younger than the least active ones. In fact, by 
inspecting the plot we can see a gradual increase in the 
equivalent width of the Li line, in going up the Figure, 
from HD 202628 to HD 174429. The spectrum of HD 
185124 was too broadened to allow some comparison with 


other stars, and was not included in the figure. The pres¬ 
ence of lithium and a high chromospheric activity are clas¬ 
sical youth indicators in late-type stars. The exceptions 
are few, and do not contradict this idea. HD 22049 (= e 
Eri) is a BY Dra variable, and possibly older than its chro¬ 
mospheric activity suggests. HD 870 has an activity level 
near that of the Vaughan-Preston gap, and could be con¬ 
sidered as an inactive star, observed during a maximum 
of activity. 

The spectra of CYKOS are presented in Fig. ||b. When 
we examine the same chromospheric activity sequence 
amongst the CYKOS, we hnd nothing similar to that 
found in normal stars. The Li line is only present in the 
spectra of HD 189931. In the others there is no trace of 
Li. The comparison can be done more properly in the bot¬ 
tom panels of this figure, where we show normal stars and 
CYKOS with the same chromospheric activity levels, and 
supposedly, the same age. 

4.3. Stellar Parameters and Models 

Stellar parameters (Teg, logg, [Fe/H]) for the program 
stars were first derived from the ubvy photometry and the 
classical relation log g* = 4.44 -|- 4 log T^/Tq -|- 0.4(Mboi — 
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Fig. 2. Spectra for active stars in the Li region (A 6707 A). Panel a presents normal young stars, while panel b 
presents CYKOS. The spectra are ordered according to the chromospheric activity, with the most active stars in the 
upper part of the figure. With few exceptions it is seen that the Li line is present in all normal active stars, and 
its equivalent width becomes smaller towards less active stars. This behaviour is totally absent from the spectra of 
CYKOS. In panels c and d we compare spectra of stars having the same activity levels and, supposedly, the same age. 
The upper spectra are for normal stars, in which there is Li. In the bottom spectra, of CYKOS, the Li line is absent. 
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4.74) -I- logM^,/M Q, usi ng the same procedure described 
in Castilho et al. (2000). 

The metallicities of the program stars were then rede¬ 
termined using curves of growth of Fe I lines where the 
updated code RENOIR by M. Spite was employed. The 
stellar parameters used, with the spectroscopic metallicity 
that we have found, are listed in Table |^. For HD 185124 
and HD 174429 a metallicity determination was not pos¬ 
sible, due to the line broadening. The error in [Fe/H] is 
0.10 dex. 

Model atmospheres employed have been interpolated 
in tabl es com puted with the MARCS code by Edvardsson 
et al. ( 1993 ). 


4.4. Li abundances 

The Li depletion must not be identical in all observed 
stars, since it depends on factors such as mass and metal¬ 
licity. Even the equivalent width of the ^Li resonance dou¬ 
blet (A 6707. SA) depends stron gly on the temperature 
of the star (Castilho et al. 200C| ). Therefore, a strong Li 
line does not always warrant a high lithium abundance. 
The stars in Fig. ^ do not have the same temperature, 
and those plots shall be taken as illustrative compar¬ 
isons. More accurate estimates of the Li age must be done 
by considering Li abundances rather than the equivalent 
widths. 

Spectrum synthesis calculations were used to fit the 
observed spectra of the stars listed in Tab. The cal¬ 
culations of synthetic spectra were carried out using a 
revised version of the code described in Barbuy (1982), 
where molecular lines of C 2 (A^H-X^H), CN red (A^H- 
X^E) and TiO 7 (A^4>-X^A) systems are taken into ac¬ 
count. The oscillator strengths adop ted are the labo rator y 
values obtaine d by Fuhr et al. (|l98^ ), Martin et al. ( 1988|) , 
Wiese et al. ( 1969 ). When these wer e not available, we 
hav e used those given by Spite et al. ( 1987 ) or Barbuy et 
al. ( 1999| ) obtained by inverse solar analysis. Solar abun- 
dances are adopted from Grevesse & Sauval ( 1999 ). The 
derived lithium abundances NH^i) of the observed stars 
are given in column 10 of Tab. 

For the determination of the lithium ages of these 
stars, we ha ve con sidered the lithium depletion diagram of 
Soderblom ( 1983 ). Fig. || presents a Li depletion diagram 
as a function of the stellar effective temperature, according 
to Soderblom (1983). The diagram is an approximation for 
the calculation of Li ages. The dashed horizontal line cor¬ 
responds to the Li abundance in the interstellar medium, 
and the curves indicate the expected depletion as a func¬ 
tion of the stellar age. Vertical arrows show the expected 
lithium depletion before the ZAMS. The vertical dotted 
lines mark the region corresponding to the Boesgaard- 
Trippico dip (cf., Boesgaard & Trippico 1986| ), where de¬ 
pletion is not linked with age, but probably with the in¬ 
ternal structure of the star. The filled squares in the top 
panel of Fig. ^ are the normal young stars of Table ||, while 
the triangles show the CYKOS. The lower panel of Fig. || 


shows the Li depletion diagram of a few CYKOS (stars) 
whose abundances are from the literature (section 5). It 
is clear that normal young stars have lithium abundances 
very similar to that of the interstellar medium. That is, in 
these stars there was no depletion. The CYKOS, on the 
other hand, have Li ages greater than 2 Gyr. 

HD 870 is the only exception amongst normal young 
stars. As mentioned before, this could be an inactive star, 
included in the sample by an error in logi?(jj^ (for in¬ 
stance, it has been observed in the H an K Ca H lines 
only once, by Henry et al. 1996| ). HD 189931, in spite of 
having N(Li) ^ 2.0, has the same Li age as HD 152391 
and HD 13445. 

We conclude that CYKOS present Li ages higher than 
the ages estimated by their chromospheric activity, al¬ 
though in better agreement with their kinematic prop¬ 
erties. 

5. Data on individual objects 


The initial suggestion by Soderblom (|1990|) that the 
CYKOS can be runaway stars seems very unrealistic. If 
these stars were young, they should present Li abundances 
typical of their youth. Their high kinematic age is a clear 
indication of their older status. An initial hypothesis that 
can be tested is that they could be chromospherically ac¬ 
tive binaries, which can be tested by searching for radial 
velocity variations. Before giving a general explanation, 
it is important to see what is known about the stars in 
Tab. 

HD 5303 = CF Tuc 

This is a GO V -I- K4 IV R S GVn-type spectroscopic 
binary (Strassmeier et al. 1993), with an orbital period of 
2.80 days. Li abundances were measured by Randich et al. 
( 1994D for both components. In Fig. ||b, both components 
are linked by a dot-dashed line. Note that this line presents 
roughly the same slope of the curves with the same Li age. 
HD 13445 = HR 637 

The index logi?(jj^ = —4.74 indicates that this star 
could be an inactive sta r. A co mparison between Ha fluxes 
(Pasquini & Pallavicini 1991 ) with other stars is ambigu¬ 
ous: HR 13445 presents fluxes similar to HD 42807 and HD 
81997 (active stars) and to HD 4308 and the Sun (inactive 
stars). Favata et al. ( |1997 ) have measured N(Li) < —0.24, 
one order of magnitude lower than our measured value. 
Queloz et al. ( 2000 ) have found a planet with 4 Jupiter 
masses at a distance of 0.11 AU from the star. 

HD 16176 = HR 756 


This is classified as F5. Balachandran (1990) measured 
N(Li) < 1.87. The depletion would be high for this spec¬ 
tral type, if due to age alone. However, the star is found 
within the Boesgaard-Trippico dip, and this low lithium 
abundance cannot be considered as indicative of old age. 
HD 20766 = Cl Ret = HR 756 

This is a G2.5 V star, visual companion of HD 2 0807, 
with an angular separation of 307". Wooley ( |l970 ) sug- 
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Fig. 3. Li depletion diagram for late-type stars, according to Soderblom ( |1983|) . Panel a: stars that we have observed; 
panel b: CYKOS with Li abundances from the literature. In this figure, the dashed horizontal line at the top indicates 
the Li abundance in the interstellar medium. Vertical arrows correspond to the Li depletion before the main sequence. 
Several curves indicate the Li depletion expected for stars of a given age. The vertical dotted lines in the left part of 
the plots correspond to the region of the Boesgaard-Trippico dip. 
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gests that they are members of the ^ Her moving group, 
which has an isochrone age of a few billion years, da Silva 
& Foy (|1987 ) measured a metallicity typical of population 
I stars and, therefore, criticized the hypothesis by Johnson 
et al. ( 1968 ) that the pair is composed of subdwarfs. There 
is no sign of radial velocity variability in either stars (da 
Silva & Foy 1987] ). 

HD 39917 = SZ Pic 

This is a RS CVn chromospherically active star 
(Strassmeier et al. 1993] ), havi ng an orbital period of 4.80 
days. However, Mason et al. (1998) have not detected a 
companion within AH < 3.0 and angular separation be¬ 
tween 0.035" and 1.08". The lithium abundance was mea¬ 


within AH < 3.0 and angular separation between 0.035'' 
and 1.08". 

HD 74385 


Dwarf K1 V. Favata et al. ( 1997 ) have measured N(Li) 
< -0.23. 

HD 88742 = HR 4013 


= 1.03. The depletion seems substantial for the tempera¬ 
ture of this star. 

HD 123651 

There is no indication of a companion within AH < 3.0 
and a ngula r separation between 0.035" and 1.08" (Mason 
et al. 199g ). 

HD 131977 

K4 V visual binary, with angular separation of 20". 
The companion is HD 1 3197 6. Its x-ray emission level is 
moderate (Wood et al. 1994|) , but there seems to be no 
doub t about its activity (R obins on, Cram & Giampapa 
1990). Duquennoy & Mayor (1988|) have found no invisible 


sured by Randich et al. (1993), and is consistent with the 
expected abundance at the ZAMS (see Fig. ||). 

HD 65721 

_ G6 V variable star, ROSAT source (Hiinsch et al. 

199S ). Mason et al. ( 1998| ) have not detected a companion 


G1 V star, ROSAT source (Hiinsch et al. 1999). Mason 
et al. (1998) have not detected a companion within AH < 
3.0 and angular separation between 0.035" and 1.08". 

HD 89995 = HR 4079 

F6 V star, ROSAT source (Hiinsch et al. 1999]) . 
The lithium abundance measured by Balachandran i 


companion for this star. 

HD 149661 = 12 Oph = V2133 Oph = HR 6171 

K2 V variable of BY Dra type (Petit 1990|) . It was 
detected by ROSAT in E UVE with moder ate i ntensity 
(Tsikoudi & Kellett 1997 ). Habing et al. ( 1996 ) report 
a Vega-like protoplanetary disk, but the presence of cir¬ 
rus during the observation has somewhat made this find¬ 
ing inconcl usive. Young, Mielbrecht & Abt ( 1987 ) and 
Tokovinin (1992) have found a constant radial velocity 
for this star, and McAlister et al. ( 1987 ) did not detect 
the presence of an unseen companion by using speckle in¬ 
terferometry. 

HD 152391 = V2292 Oph 


N(Li) = 2.38) is high compared to some stars, but small 
for the temperature of this star (6280 K). The star is lo¬ 
cated within the Boesgaard-Trippico dip, where depletion 
is uncorrelated with age. 

HD 103431 

It is a dG7 star, visual companion of HD 103432, an¬ 
gular separation of 73.2". Gonstant radial veloci ty du ring 
a time span of 2499 days (Duquennoy & Mayor 1991 ). 

HD 106516 = HR 4657 

F5 V spectr oscop ic binary with period of 853.2 days 
(Latham et al. 1992]). Li thium abundance was measured 
by Lambert et al. ( 1991 ), N(Li) < 1.32. The star is also 
located with in the Boesgaard-Trippico dip. Edvardsson 
et al. ( 1993 ) calculate an i sochrone age of 5.37 Gyr. 


Fuhrmann & Bernkopf ( |l999|) suggest that this star is a 
field blue straggler, having a chemical composition and 
kinematics typical of thick disk stars, in spite of having 
an age typical of thin disk stars. 

HD 120237 = HR 5189 

It is classified as G3 IV-V. There is no indication of a 
companion within AH < 3.0 and ang ular separation be¬ 
tween 0.035" and 1.08" (Mason et al. 199^ ). The lithium 
abundance was calculated by Randich et al. ( 1999|) , N(Li) 


G8 V variable star of BY Dra type (Petit 1990 ). 
Detected by ROSA T in E UVE with moderate intensity 
(Tsikoudi & Kellett 1997 ). Gonstant radial velocit y dur - 
ing a timespan of 3387 days (Duquennoy & Mayor 1991 ). 

HD 165401 

It is a GO V, relativaly metal-poor star ([Ee/H] ~ 
—0.50 dex). The index logi?(jj^ that we h ave used refers 
to a sole observation (Duncan et al. 1991 ). We have con¬ 
sidered the possibility that this star is inactive, but its 
emission in Ha (Herbig 1985) seems consistent with its 
logi?(jj^. There is no unseen companion with a magni¬ 
tude diference lower than 2.5 mag and angular separation 
greater th an 3 AU, according to sp eckle interferometry 
(Lu et al. 1987 ). Abt & Levy ( 1969 ) have found a radial 
velocity variability of ±5 km/s, but recent investigations 
do not confirm this result (Duquennoy & Mayor 1991; Abt 
& Willmarth 1987). Guriously, radial velocities measured 
for this star during the sixties and seventies yield values 
homogeneously around Vr ~ —114 km/s, while all recent 
studies yield values around Vr « —120 km/s. 

HD 196850 

G2 V star, with constant radial velocity during a time 
span of 3995 days (Duquennoy & Mayor 1991). 

BD -f 15 3364 

Solar-metallicity GO V star. Duquennoy & Mayor 


(1991) have found a constant radial velocity during a time 
span of 2264 days. Garney ( |1983 ) found no photometric 
variability characteristic of an unseen companion. 
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6. The nature of CYKOS 

6.1. Coalescence of close binaries 

The majority of objects considered in the previous section 
are undoubtedly active, have lithium abundances lower 
than that of stars with 2 Gyr of age, and some show no 
indication of radial velocity variations. They are thus sin¬ 
gle objects, most probably old. 

As we have seen, some of them are chromospherically 
active binaries. Their chromospheric activity results from 
the synchronization of the orbital with the rotational mo¬ 
tion. This is the reason why they have a low chromospheric 
age, but a higher age from the point of view of their kine¬ 
matics. Nevertheless, in the case of single stars, there is 
no known mechanism that could store angular momentum 
to be used later by the star. The anomalously low lithium 
abundance, together with the high velocity components, 
are hardly interpreted as a consequence of anything but 
an old age. Even the BY Dra variables in the sample must 
be single stars, and not chromospherically active binaries, 
since there is no indication of variability in Vr for them. 
Amongst BY Dra stars are found binaries as well as single 
stars (Eker 1992). Soderblom ( 1990| ) has estimated a kine¬ 
matic age of 1-2 Gyr for them, contrary to the expectation 
that these stars could keep high activity levels at advanced 
ages. The low kinematic age must be understood as an ev¬ 
idence that the BY Dra-type variability is not exclusive 
of chromospherically active stars, but could be generated 
by the intensity of the magnetic activity itself amongst 
low-mass stars, binaries or not. 


used by Fuhrmann & Bernkopf (1999). 

The coalescence scenario was already considered as a 
classic explanation for the formation of blue stragglers. 
Two works (van’t Veer & Maceroni 1989; St§pieh 1995) 


Poveda et al. (1996a ,b) also have found several CYKOS 
amongst UV Ceti stars, which are known as young low- 
mass stars. The authors suggest that these objects could 
be red stragglers, a low-mass analogous of the well-known 
blue stragglers. According to this hypothesis, the red 
stragglers would be produced by the coalescence of two 
low-mass stars (with about 0.5 Mq each) originally in a 
short-period binary pair. According to them, these ob¬ 
jects would be stragglers in a velocity diagram, compared 
to other stars. However, being more rapid, they would not 
be stragglers in the sense they are in the HR diagram. 
Thus, the name ‘held blue stragglers’ includes not only 
the idea about their origin, but also their present location 
out of a cluster. Note that the same denomination was 


that investigate the coalescence of short-period binaries 
into a single star predict the formation of low-mass blue 
stragglers. St§pieh (1995) has even shown that the coales¬ 
cence is more easily attained for low-mass binaries (each 
having around 0.6 Mq, forming a star with 1.2 Mq) than 
for more massive stars that originate the classic blue strag¬ 
glers in open clusters. The formation of a low-mass blue 
straggler could be achieved within 2.5 Gyr, for systems 
having an initial orbital period of 2 days. 


The properties of a supposed low-mass held blue strag¬ 
gler would be tightly similar to that of some GYKOS, as 
we will see in what follows. 

In short-period binaries, we expect the occurrence of 
synchronization between the orbital and rotational peri¬ 
ods. For low-mass stars, the magnetic activity is strong, 
and increases the angular momentum loss. When both or¬ 
bital and rotational periods are synchronized, the rota¬ 
tional angular momentum loss occurs at the expenses of 
the orbital angular momentum. As a result, the period de¬ 
creases, the components rotate more rapidly, and become 
closer to each other, eventually becoming contact binaries, 
as those of W UMa-ty pe. 

Rasio & Shapiro ( 1995| ) simulate systems like these, 
using the technique of smooth particle hydrodynamics. 
The authors show that once the contact is achieved, these 
systems are dynamically unstable and can rapidly coa¬ 
lesce into a single object, having a high rotation rate. The 
events related to this coalescence can produce an extense 
outer envelope that would make the star appear like a 
pre-main sequence star. According to the authors, the co¬ 
alescence can occur in a time scale of a few hours, once the 
dynamical instability is set. The envelope is kept gravita¬ 
tionally bound to the star, which eventually contracts to¬ 
wards thermal equilibrium. Mass loss, in this event, would 
be minimal. 

In a coalescence of two low-mass stars (0.5 Mq each), 
the resulting star must present a mass similar to that of 
the Sun and a high rotation rate. This rotation rate, to¬ 
gether with the convection in the outer stellar atmosphere, 
would produce a copious chromospheric activity, similar to 
the one found in very young stars. In the case of low-mass 
stars that have not ignited hydrogen in their cores, the just 
formed single star would be similar in many respects to a 
young star, positioning in the zero-age main sequence, like 
the blue stragglers. However, this star would inherit the 
same velocity components of the binary pair from which it 
was formed. Thus, due to the time before the coalescence, 
the velocity components are not similar to that of a young 
star. We would have a star almost in everything young, 
but kinematically old. The lithium abundance is also one 
of the few tracks that can show the real nature of these 
objects. In spite of not burning hydrogen considerably, 
stars having around 0.5 Mq are highly convective, and 
Li burning is very efficient in them. Blue Stragglers like 
these should present small or no Li abundance (Pritchet 
& Glaspey 1991 ; Glaspey et al. 1994 ), since they would be 
formed by older objects. 

A criticism that could be made is that, if CYKOS 
are field blue stragglers, formed during the coalescence of 
low-mass short-period binaries, why do their spectra not 
present very broad lines, as it is expected in stars with 
high rotation rates? The same problem occurs for blue 
stragglers in open and globular clusters, that do not rotate 
more rapidly than young normal stars. Leonard & Livio 
( 1995 ) suggest that the majority of the angular momen¬ 
tum is stored in the extense disk that is formed around just 
formed blue stragglers. The central object would expand. 
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due to the thermal energy resulting from the coalescence, 
and contracts toward thermal equilibrium, in a time scale 
lower than 10^ years, similar to that of pre-main sequence 
stars. Thus, we do not find CYKOS with broader lines 
simply because we do not expect that they all have been 
formed within the last 0.5 Gyr. 

Seven amongst our stars seem to fit well within these 
criteria: HD 20766, HD 106516, HD 131977, HD 149661, 
HD 152391, HD 165401, BD -1-15 3364. Other stars are 
suspect, but the information is insufficient to characterize 
them as low-mass blue stragglers. One of the stars (HD 
189931) seems to be a real young single object with high 
velocity components (see also Tab. |^). Nevertheless, there 
is little information in the literature that could test this 
hypothesis. The other stars need more data to investigate 
whether they are field blue stragglers, chromospherically 
active stars or runaway stars. 

6.2. Coalescence Rate 


that T^pm can be approximated by a function Tp, which 
depends only on the initial period, and that is given by 
the time for contact for the system O.QMq + 0.6Mq. Note 
that being rp a maximum time scale, our estimates will 
be a little underestimated, since the total number of sys¬ 
tems that have achieved contact rp after the formation 
of the binary will be greater that the calculated number, 
due to the number of binaries with lower total masses that 
achieve contact more rapidly. Also, due to the initial mass 
function, the number of systems with lower total masses 
must be higher. Taken these into account, Eq. reduces 
to Nb(t) ^ cF{t), where 

/.1.2 

c = b (j){m) dm (3) 

Jo.8 

and 

poo 

P{t)=i Tp)d\ogP, (4) 

Jo.3 


The formation rate of these objects can be calculated from 
considerations about the initial mass function, star for¬ 
mation rate, initial period distribution in binaries and the 
time scale for contact. The time needed for the coales¬ 
cence, once contact is achieved, is very smal l com pared to 
the time scale for contact (Rasio & Shapiro 1995 ), and we 
will consider it negligible. 

It can be shown that the total number of stars, having 
masses between mi and m 2 , that becomes contact binaries 
in t is 


rO.5 


Nb{t) = b 


Pin, logP, m, t) d/xdlogPdm(l) 


/0.3 ^0 


'P{p,,\ogP,m,t) = f{^i)(p{\ogP)il}{t - (2) 

where r^pm is the time scale for contact for a binary, hav¬ 
ing initial period in the range (log P, log P + AlogP), 
secondary-to-total mass ratio between {p, ^ -|- A/r) and 
total mass between (m,m + Am), (^(logP) is the initial 
period distribution, /(/r) is the distribution of mass ra¬ 
tio, '0(t) is the star formation rate, and (f>{m) is the initial 
mass function. 

For the computation of the equation above, we will 
consider i^(logP) « 0.018 -I- 0.027 log P, with a cutoff 
for logP < 0.3, which approximates fairly well the ini¬ 
tial period distribution of Duquennoy & Mayor (1991), 
for the region 0.3 < logP < 1.0, where P is in days. We 
have assumed f{fj.) = 24/r^ (Matteucci & Greggio 1986D , 
mi = O.SMq and m 2 = I. 2 M 0 , which correspond to pairs 
having total masses equal to the mass range for G dwarfs. 

The function r^pm is very complicated. St§pieh (1995) 
published calculations for only three binary configura¬ 
tions: IM 0 -1- IM 0 , IM 0 -|- O.65M0 and 0.6Mq + O.GMq. 
The last of these has total mass equal to m 2 , but calcula¬ 
tions for lower total masses are not published, neither for 
different /r. However, St§pieh says that configurations for 
lower total masses achieve contact in a time scale lower 
than that for the system 0.6Mq + 0.6Mq. We consider 


since the integral in /i is unity. 

For the calculation of rp, we have considered Fig. 2 of 
St§pieh ( 1995| ). The time scales for contact were fitted by 
Tp = 0.8e^log^ have assumed a constant star forma¬ 

tion rate, in spite of the evidences for its non-constancy 
(Rocha-Pinto et al. 2000b), since we are only interested in 
the magnitude of the coalescence rate. 

In Fig. IJa, the initial period distribution is shown. The 
time scale for contact is shown in panel b. The number of 
field blue stragglers already formed was calculated for a 
sphere of radius 30 pc around the Sun, which correspond 
approximately to the volume within which our sample is 
nearly complete. This rate is shown in Fig. ^c. Integrating 
the data in this panel, we have a total number of 28 blue 
stragglers already formed in this sphere. 

The calculations are very rough, as can be seen from 
the approximations considered. However, the number of 
objects formed in the sphere is similar in order of magni¬ 
tude to the number of field blue stragglers in our sample 
(including suspect objects). This reinforces our hypothesis 
for the nature of these objects. 

From these results, we believe to hav e foun d popula¬ 
tion I field blue stragglers. Poveda et al. ( |l996a| ,b) arise a 
hypothesis that cannot be tested with the same probabil¬ 
ity level as ours, since they do not have an indication of 
old age for their stars, besides their kinematics. 

The traditional method used for finding a blue strag¬ 
gler, based on its position in the HR diagram, is not pos¬ 
sible to be used for a population I field blue straggler. 
The search for objects with strong chromospheric activ¬ 
ity, or little Li, or with high velocity, also do not allow 
their identification, since normal stars present these prop¬ 
erties. Only by the intersection of several properties we 
have found objects which apparently can be best explained 
by this scenario. 

It is worth mentioning an independent research on 
ultra-lithium -defici ent halo stars recently published by 
Ryan et al. ( 2001 ), which conclude that such stars can 
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Fig. 4. Calculations for the formation rate of population I field blue stragglers having masses between 0.8-1.2 Mq. 
Panel a: initial period distribution, adapted from Duquennoy & Mayor ( |l99l| ); panel b: time scale for contact; panel 
c: number of field blue stragglers, with mass in the considered range, formed within a sphere of 30 pc radius centred 
at the Sun. 


have the same origin of the blue stragglers, being their 
lower mass the only significant difference between them. 
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